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ABSTRACT 

We present spectroscopic observations in the rest-frame optical and near- to mid-infrared wave- 
lengths of four gravitationally lensed infrared (IR) luminous star-forming galaxies at redshift 1 < z < 
3 from the LUCIFER instrument on the Large Binocular Telescope and the Infrared Spectrograph on 
Spitzer. The sample was selected to represent pure, actively star- forming systems, absent of active 
galactic nuclei. The large lensing magnifications result in high signal-to-noise spectra that can probe 
faint IR recombination lines, including Paa and Bra at high redshift s. The sample was augmented 
by three lensed galaxies with similar suites of unpublished data and observations from the literature, 
resulting in the final sample of seven galaxies. We use the IR recombination lines in conjunction 
with Ha observations to probe the extinction, Av, of these systems, as well as testing star formation 
rate (SFR) indicators against the SFR measured by fitting spectral energy distributions to far-IR 
photometry. Our galaxies occupy a range of Av from ~ to 5.9 mag, larger than previously known 
for a similar range of IR luminosities at these redshifts. Thus, estimates of SFR even at z ~ 2 must 
take careful count of extinction in the most IR luminous galaxies. We also measure extinction by 
comparing SFR estimates from optical emission lines with those from far-IR measurements. The com- 
parison of results from these two independent methods indicates a large variety of dust distribution 
scenarios at 1 < z < 3. Without correcting for dust extinction, the Ha SFR indicator underestimates 
the SFR; the size of the necessary correction depends on the IR luminosity and dust distribution 
scenario. Individual SFR estimates based on the 6.2 /im PAH emission line luminosity do not show a 
systematic discrepancy with extinction, although a considerable, ~ 0.2 dex scatter is observed. 

Subject headings: cosmology: observation — galaxies: evolution — galaxies: high-redshift — infrared: 
galaxies 



1. INTRODUCTION 

The evolution of the star formation rate (SFR) of 
galaxies is a central topic to the study of galaxy evo- 
lution. It is generally agreed that the SFR density of 
the Universe has declined by a n order of magnitude 
since z ^ 1 to the p resent fe .g.. lLe"Floc'h et al l 120051 
lR.uiopakarn et al.l[201 0: Magn elh et al.ll2011[ ). The exact 
epoch of the peak of the SFR history is not known pre- 
cisely, although it app ears to be constrained to be within 
1 < z < 3 (e.g.,|P erez- Gonzalez et al. 2005; Re ddv et al.l 
120081 : iRodighiero ~h 2010: Magnelh et al 20m 

At this redshift range, the primary SFR indicators are 
based on infrared (IR), optical, and extinction-corrected 
ultraviolet (UV) observations; the resulting SFR esti- 
mates commonly disagree with each other by more than 
a factor of two (e.g., Reddy et all l2008f ). The major- 
ity of star formation at these redshifts is known to oc- 
cur in optically extinc ted star-forming regions in IR lu- 
minous g:alaxies (p-R-. Le"Floc'h et al.ll2005l iDole et aTl 
120061 : iBerta et al.H2011f ). and the uncertainties resulting 
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from the extinction undermine our ability to understand 
the SFR history during this era. A better understanding 
of the SFR history will have important implications on 
the cosmic stellar mass build-up and metal production. 
For example, there is a significant discrepancy between 
the expected metal abundance derived fro m the SFR his - 
tory and the observed abundance (Bou chJet al.l 120071 ). 
with some studies indicating the difference to be nearly 
an order of magnitude, and as a result placing the peak 
of the SFR history as far back as z ~ 4 (jKobavashi et al.l 

The situation necessitates more efforts into exploring 
the nature of optical extinction in high redshift galaxies 
using unbiased measures as well as studying the effect of 
extinction on various SFR indicators. 

Extinction measurements based on optical emission 
lines could be biased in highly obscured environments 
because the indicator only probes the outer layer of 
the star- forming regions, where the extinction is rel- 
atively low. This effect is observed locally (e.g., 
lAlonso-Herrero et al.ll2006l ). but only because observa- 
tions are available for longer wavelength IR recombina- 
tion lines (e.g. Paa and Bra) that are less affected by 
extinction. In this work, we extend this technique out 
to z = 3 by comparing the strength of the Ha line with 
those of Paa and Bra. Since the latter are in wavelength 
regions with 10 — 20 times less extinction than Ha, they 
provide a measurement of extinction through the entire 
star-forming region. With an unbiased estimate of ex- 
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tinction, we can address its effects on SFR indicators in 
this critical redshift range. 

In this paper, we study a wide range of star formation 
diagnostics in seven gravitationahy lensed star-forming 
galaxies at 1 < z < 3. We observed four of these galaxies 
spectroscopically with the LUCIFERl near-IR imaging 
spectrograph on the Large Binocular Telescope (LBT) 
to measure the Ha line flux in the near-IR, and with the 
Infrared Spectrograph (IRS, 5 — 38 /im) on the Spitzer 
Space Telescope to observe the wavelength regions cov- 
ering Paa and/or Bra lines as well as aromatic emis- 
sion lines and emission complexes (commonly attributed 
to, and hereafter, polycyclic aromatic hydrocarbons or 
PAH). The sample comprises Abeh 2218b, Abeh 2667a, 
Abeh 2218a, and Abeh 1835a at redshift 1.03, 1.03, 2.52, 
and 2.57, respectively (magnification 6x - 27x). The 
IRS observations at longer wavelengths are further used 
to compare SFR estimates from the 6.2 /im PAH fea- 
ture. They were carried out under Spitzer program ID 
82, 30775, 50586 (PI G. Rieke); and 40443, and 50372 (PI 
C. Papovich). Our sample is augmented by three galax- 
ies with similar suites of observations from unpublished 
data and the literature: SDSSJ120601+5142 (hereafter 
the Clone), the Lyman break galaxy LBG MS 1512-cB58 
(hereafter cB58), and the 8 O'clock arc, at redshift 2.00, 
2.73, and 2.73, respectively. Our galaxies were also ob- 
served with the MIPS instrument (Rieke et al.""2004) on 
Spitzer at 24 and 70 /im by Rigby et al. (2008), allow- 
ing us to combine the MIPS data with far-IR photometry 
from the literature to estimate SFR via I/(TIR) obtained 
by fitting the spectral energy distribution (SED) across 
the peak of the far-IR emission. In addition, we used 
the 2 4 /im monochromatic indicator (jRujopakarn et aP 
[20121 ) to estimate L(TIR). 

This paper is organized as follow. We discuss the LBT 
and Spi^^^er observations and data reduction in Section [21 
extinction, metallicity, L(TIR), and SFR measurements 
in Section [3l compare SFR indicators in Section [H and 
make concluding remarks in Section [5] We assume a 
ACDM cosmology with = 0.3, I^a = 0.7, and Hq = 
70 km s~^Mpc~^ throughout this paper. 

2. OBSERVATIONS AND DATA REDUCTION 

Here we describe the selection of our sample 
(Section 12. ip : the spectroscopic observations with 
LBT/LUCIFER along w ith t he near-IR data reduction 
(Section 12.21 and Section 12. 3p and with Spitzer (Section 
12. 4p : data compilation from the literature (Section [23]) : 
additional photometric observations with Spitzer MIPS 
and the respective data reduction procedure (Section l2.6p 
that results in the final data set for our sample of galax- 
ies. 



2.1. Sample Selection 

To measure extinction and study its implication on 
the SFR indicators, we need a sample of star-forming 
galaxies that we understand very well. Gravitationahy 
lensed galaxies are outstanding candidates because their 
amplified fluxes and images enable high signal-to-noise 
(S/N) spectroscopy. The candidates for our new obser- 
vations were drawn from the pool of objects studied by 




Fig. 1. — Images of the lensed galaxies in our sample from the 
Hubble Space Telescope Advanced Camera for Survey (ACS) with 
the positions of our LBT/LUCIFER slit for NIR spectroscopy 
shown in red. Clockwise from top-left are Abell 2218b, Abell 2218a, 
Abell 1835a, and Abell 2667a. The ACS filter for Abell 2218a and 
Abell 2218b is F625W; Abell 1835 and Abell 2667 images were 
taken in F850LP. Each image is 37^^ wide (identical scale in all four 
panels), north is up, east is left. For Abell 2218a and Abell 26 67a, 
each lensed component is labeled following IKneib et al.l (|2004l ) and 
ISand etldl (|2005 V respectively 

iRigbv et all (j2QQ8f ). Briefly, the lRigbv et~all (2008) selec- 
tion requires that objects have (1) Spitzer MIPS 24 /im 
flux above 0.4 mJy to be observed efficiently with IRS; 
(2) morphologies in the optical iJiST imaging that exclude 
members of the lensing cluster; (3) spectroscopic or prob- 
able photometric redshift a bove 1; (4) lensing magnifica- 
tion above 3x. From the Rigb v et aTl (jlOOS) sample, we 
further require that candidates are absent of AGN ac- 
tivity based on X-ray luminosity below 10^^ erg/s/cm^ 
and the lack of mid-IR power-law SEDs, indicating that 
their 24 jam fluxes are not dominated by AGN emiss ion. 
Apart from the objects selected from the Rigby et al.l 
(|2008l ) sample, we include three additional objects that 
match the selection criteria; cB58, the Clone, and the 8 
O'clock arc. We will discuss these objects in Section [231 
The galaxies in our sample are located behind galaxy 
clusters that are well-modeled for mass distribut ion. 
The mass model of the Abell 2218 cluster (jKneib et al.l 
Il996[ ) indic ates a magn i ficatio n value for Abell 2218a 
of 22 ± 2 (|Kneib et al.l [200l . For Abeh 2218b and 
Abell 2667a, magnification estimat es of 6 .1 and 17, re- 
spectively, are modeled by Rigbv ^ al.l ([2008). For 
Abell 1835a , we adopt the magnification estimate from 
iSmail et al.l (j2005l 3.5 ± 0.5. The sample is listed in 
Tabled 

2.2. LBT/LUCIFER Near-Infrared Spectroscopy 

Our near-IR spectroscopy was obtained in three cam- 
paigns during 2010 October and 20 11 May using the LU - 
CIFERl instrument on the LBT (jMandel et al.ll2007[ ). 
The common aspects of all these observations are that we 
use the N1.8 camera with plate scale of 0'.'25 per pixel and 
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TABLE 1 

The Sample of High-Redshift Strongly-Lensed Star-Forming Galaxies 



Source 


z 




R.A. 


Decl. 


Lensing 


Reference 






(J2000) 


(J2000) 


Magnification 






Abell 2218b 


1.034 


16 


35 55.16 


+66 11 50.8 


6.1 


1, 


2 


Abell 2667a 


1.035 


23 


51 40.00 


26 04 52.0 


17 


1, 


2 


The Clone 


2.003 


12 


06 01.80 


+51 42 30.7 


27+1 


3, 


5 


Abell 2218a 


2.520 


16 


35 54.18 


+66 12 24.8 


22 + 2 


1, 


2 


Abell 1835a 


2.566 


14 


01 04.96 


+02 52 24.8 


3.5 + 0.5 


1, 2 


, 6 


cB58 


2.729 


15 


14 22.29 


+36 36 25.7 


- 30 


4, 


7 


8 O'clock 


2.731 


00 


22 40.97 


+14 31 14.0 


- 8 


8, 


9 



Note. — Refe rences: (1) redshifts from our LE T Ha spectroscopy ; (2) 
IRigbv et al." ("2008"); (3) "Hainline et al. (200 a): {^^ iTeplitz et al.l (12000): f5) 
[Hn et al. (2009); (6) Small et al. (2005); (7) Seitz et al.l (|l998l ): (8) lAllam et al.l 
(|20071 ): iFinkelstein et al.. (,2009i ) 



with the 210 hnes/mm high-resolution grating, resulting 
in spectral resolution of 5.0 A. A-stars were observed for 
spectrophotometric calibration for each target (details 
on spectral type below) at similar airmass as the science 
exposure, either immediately before or after the science 
exposures. Positions of the LUCIFER slit for each of our 
observations are shown in Figure [TJ 

Abell 2667a was observed on 2010 November 3 in the 
longslit spectroscopy mode using a slit of V% by 3.9 ar- 
cmin under variable seeing of 1'.'2 — 2'.' 5. The slit was 
rotated to P.A. = 320° to place it along the length of 
the lensed arc. The total integration time was 12 x 300 
s with the telescope nodded 50'' along the slit. The un- 
usually long nodding was required by the ~ 25'' length 
of the arc. For calibration, we observed HD 223466, an 
A3V calibration star with V = 6.42 mag, through the 
same setup using a total integration time of 2 x 20 s in 
two nodding positions. 

Abeh 2218b was observed on 2011 March 10 in the 
longslit spectroscopy mode using a slit of I'.'O by 3.9 ar- 
cmin under 0'.'4 — 0'.'5 seeing. The slit was rotated to P.A. 
= 42.5° to place it along the lensed arc, which is ~ 20" 
long. The arc is 1'.'9 wide in the widest region; thus we 
took the spectra at the central position and another po- 
sition shifted by 0'.'9, with the O'.'l overlapping allowing 
for possible pointing errors to prevent a gap in the slit 
mapping. We will refer to these positions as the "Cen- 
tral" and the "Central+0'.'9" positions, respectively. The 
total integration time was 6 x 300 s at the Central posi- 
tion and 5 x300 s at the Central+0'.'9 position, totalling 
55 minutes, with the telescope nodded 40" along the slit 
for both. An AOV calibration star, HD 145454, with V 
= 5.43 mag was observed using a total integration time 
of 2 X 20 s in two nodding positions. 

Abeh 2218a was observed on 2011 March 11 in the 
longslit spectroscopy mode using a slit of I'.'O by 3.9 ar- 
cmin under 0'.'9 — I'.'l seeing. The slit was rotated to P.A. 
= 18.9° to place it along the lensed arc. The total inte- 
gration time was 10 x 300 s with the telescope nodded 
20" along the slit. We observed the same calibration star 
as for Abeh 2218b. 

Abell 1835a was observed on 2011 May 7 in the longslit 
spectroscopy mode using a slit of 1'.'5 by 3.9 arcmin un- 
der I'.'O -I'.'l seeing. The slit was rotated to P.A. = 15.0° 
to place it along the length of Abell 1835a. The total in- 
tegration time was 12 x 300 s with the telescope nodded 
20" along the slit. An A2V calibration star, HD 122365, 



with V = 5.98 mag was observed using a total integration 
time of 2 X 20 s in two nodding positions. 

We took Ar lamp exposures for wavelength calibration 
along with dark and fiat exposures during daytime. 



2.3. LBT/LUCIFER Data Reduction 

Our near-IR spectral reduction has four steps: sky sub- 
traction, spectral extraction, flux calibration, and line 
flux measurements. We first use a modified version of an 
IDL longslit reduc tion package wr itten by G. D. Becker 
for NIRSPEC (Be cker et ani2009f ) t o perform s k y sub - 
traction following the prescription of iBian eTall ^OW). 
The sky-subtracted 2-dimensional spectra are corrected 
for the 2D dispersion distortions, wavelength calibrated, 
and extracted using IRAF. We note that the flexure 
compensation mechanism of LUCIFERl was not avail- 
able during our campaigns, but clean sky subtraction 
was achieved with the reduction procedure. We com- 
bine spectra for each object by averaging values within 
each wavelength element with sigma clipping at 2.5cr 
and then flux calib rate using the xtellcor^eneral tool 
(|Vacca et al.l 120031 ) and the A- star calibrators described 
in the previous section. Line fluxes for our objects are 
measured by integrating the line and subtracting the 
background estimated from the mean continuum in the 
spectral range. We estimate the line flux uncertainties by 
a Monte Carlo simulation done by repeatedly (n = 10^) 
simulating the object's spectrum using random values 
drawn from Gaussian distribution centered at the ob- 
served flux and a equal to the observed uncertainties 
and repeat the line measurement. We take the resulting 
l-<j distribution of the simulated line flux values to be 
the uncertainties for the line flux. 

2.4. Spitzer IRS Mid-Infrared Spectroscopy and Data 
Reduction 

Unpublished Spitzer IRS spectroscopy in this work is 
from: (1) Spitzer program 50586 (PI G. Rieke) to take 
deep (integration times ~15 ks) first order Short-Low 
spectra (SLl; 7.4 — 14.5 jam) in the Bra wavelength re- 
gions of Abell 2667a and Abeh 2218b; (2) Spitzer pro- 
gram 40443 (PI C. Papovich) to take deep (integration 
time 7—14 ks) second order Short-Low spectra (SL2; 
5.2 — 7.7 /im) in the Paa wavelength regions of Abell 
2218a and Abell 1835a. For each of the four objects. 
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TABLE 2 

Near- and Mid-Infrared Spectroscopic Observations 



Source 


IRS Exposure Time (ks) 


IRS Program ID 


Near-IR Spectroscopy 




SL2 


SLl 


LL2 


LLl 






Abell 2218b 




15.1 


3.6^ 


3.6^ 


30775, 50586 


This Work 


Abell 2667a 




16.1 


1.9^ 


1.9^ 


30775, 50586 


This Work 


The Clone 


8.6 




21.6^ 


20.9^ 


40430, 50372 


Hainline et al. (2009) 


Abell 2218a 


14.4^ 






3.6^ 


82, 40443 


This Work 


Abell 1835a 


7.7 






3.6^ 


82, 40443 


This Work 


cB58 


6.8^ 


6.8^ 


14.6^ 


34.1^ 


30832 


TeDlitz et al. f2000) 


8 O'clock 


14.4 






22.1 


40443 


Finkelstein et al. (2009) 



Note. — Spit zer IRS spectra previousl y pu blished in ( a) Ri gbv et al. (2008); (b) 
IFadelv et all (f2010l ): (c) IPapovich eFall ((20091 ): (d) ISiana et al.1 (f2008i )" 



the IRS first order Long-Low spectra (LLl; 19.5 — 38.0 
/im) were from GTO programs 82 (Abell 2218a and Abell 
1835a) and 30775 (Abell 2218b and Abell 2667a), PI G. 
Rieke. The IRS second order Long-Low spectra (LL2; 
14.0 - 21.3 /im) are available for Abeh 2218b and Abeh 
2667a from program 30775, and Abel l 1835a from pro - 
gram 82. These data are published in iRigbv et al.l ()2008l ). 
The details of IRS observations for our sample are sum- 
marized in Table [21 

We have reduced all archival spectroscopic data for ob- 
jects in our sample using the most updated version of the 
processing software from L evel- 1 data (BCD) fo r IRS by 
following the prescription of lTeplitz et al.l ()2007l ) . Briefly, 
first we use IRAF to fit the background slope with time 
for each row and subtract that fitted value row by row 
to remove latent charge. Second, we interpolate over 
bad and/or hot ("rogue") pixels using the IDL routine 
IRSCLEAN (version 2.1) provided by the SSC based on 
the known hot pixel mask for the corresponding cam- 
paign, plus manual identification. Third, we subtract 
sky using a sky image constructed from a median of other 
map positions. Then we co-add each map position into 
a 2-dimensional (2D) spectrum and use the SPICE soft- 
ware (version 2.4) provided by the SSC to extract a 1- 
dimensional (ID) spectrum using the optimal extraction 
template for point sources for each map position. The 
ID spectra are combined with sigma clipping at 2.5cr into 
the final ID spectra shown in Figures [2] and [3l 



line measurements are determined with this iterative fit- 
ting procedure. This method of continuum estimation al- 
lows construction of a mid-IR spectral template that best 
matches each galaxy in a non-parametric manner. In 
other words, it provides continuum estimates in the Paa 
and Bra wavelength regions that are not only constrained 
by photometric observations (e.g., from S pit zer/IR AC), 
but also by PAH-region information. These continuum 
estimates are shown for each object in Figures O [3l and 

SI 

Bra and Paa lines are measured by fitting Gaussian 
profiles with widths fixed at the resolution of the IRS 
module covering the line, the central wavelengths fixed 
at the expected line wavelengths based on the optical 
spectroscopic redshift, and the Gaussian base given by 
the continuum under the expected wavelength range of 
the line; the line peak is the only free parameter to fit. 
The line flux uncertainties are estimated by a procedure 
similar to that of the Ha line flux uncertainties: (1) a se- 
ries of spectra (n = 10^) were generated randomly based 
on the actual spectra and the a uncertainties; (2) line 
flux measurements are carried out on them; and (3) we 
take the l-cr value of the simulated flux distribution to 
be the uncertainties for the line flux. If the formal fit 
to the line yields negative line flux due to non-detection 
(occurring in two cases, cB58 and The Clone), we adopt 
the l-cr value as an upper limit to the line flux. The 
measurements are reported in Table [3l 



We use the PAHFIT package ( Smith et al.l I2007D to 
measure PAH feature fluxes. PAHFIT uses minimiza- 
tion to simultaneously fit the PAH and nebular emission 
lines, stellar and dust continua, as well as the silicate 
absorption feature. We report the PAH flux measure- 
ments in Table m To define the dust and stellar contin- 
uum for the measurement of Paa and/or Bra line fluxes, 
we iteratively fit the IRS spectra with PAHFIT and in- 
terpolate over points with fit residuals greater than 3cr. 
This process is repeated until there are no remaining out- 
liers (usually within the 3^^ iteration) and the final fit, 
which is effectively a noise- free template representing the 
object's spectrum, is adopted as a continuum estimate. 
The actual line flux measurements for both the PAH 
emissions and recombination lines are done on the origi- 
nal spectra, i.e. only the continua for IR recombination 



2.5. Archival and Literature Data 

We augmented our sample with three additional ob- 
jects: The Clone, cB58, and the 8 O'clock arc. These 
objects have suites of data similar to our sample: rest- 
frame optical spectroscopy covering Ra and deep IRS 
spectros copy cover ing Paa (Table [T]). 

cB58 (jYee et al.l 1 1996). at redshift z = 2.729, is the 
first lensed LBG known with high magnification (~30x); 
it has been studied extensively from UV to millime- 
ter wavelengths (Siana et al. 2008, and the references 
the rein). Rest-frame o p tical spectroscopy was obtained 
bv lTeplitz et al.l (j2000f ). ISiana et al.l (j2008f ) present a fuh 
suite of IRS observation (SL2, SLl, LL2, and LLl) for 
the galaxy. We re-reduced these IRS data using our pro- 
cedure above and found the result to be consistent with 
the original reduction published in .Siana et al.. (20081 ): 
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Fig. 2. Spitzer/IRS mid-IR and LBT/LUCIFER near-IR spectra of Abell 2218b. Panel (a) is the magnified region of IRS SLl spectrum 
at the Bra wavelength. The best Gaussian fit to the Bra emission line is the thick black line. The color-coded Gaussians at the Bra line 
wavelength illustrate the expected Bro; line fiux for each extinction scenario from Av of to 8.0 mag given the observed Ha fiux. Panel (b) 
shows the Ha emission line from LBT/LUCIFER. The light and dark gray spectra are extracted from the two parallel slit positions shown 
in Figure ^ A double-peak line profile is clearly visible. The continuum level for each aperture is shown as dashed lines. Panel (c) presents 
the overall IRS spectra with the green, orange, and red lines representing spectrum from the IRS SLl, LL2, and LLl modules, respectively. 



TABLE 3 

Measured Fluxes 



Source 


/(24 fim) 


/(70 fim) 


/(Ha) 


/([Nil] 6583) 


/(Paa) 


/(Bra) 


/[Ne 


II 12.8] 


/[Ne HI 15.5] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




(7) 


(8) 


Abell 2218b 


1.67 


7.4 ± 1.5 


21.2 ± 0.6 


5.2 ± 0.5 




4.3 ± 0.9 


1.8 


± 0.2 


0.9 ± 0.2 


Abell 2667a 


4.52 


19.4 ± 3.9 


78.7 ± 1.5 


31.1 ± 3.2 




0.4 ± 1.0 


2.8 


± 0.4 


0.9 ± 0.3 


The Clone 


0.88 


8.8 ± 1.8 


20.2 ± 0.5^ 


3.9 ± 0.2^ 


< 1.3 










Abell 2218a 


1.16 


< 7 


6.2 ± 0.3 


1.5 ± 0.3 


5.7 ± 1.2 










Abell 1835a 


0.99 


< 13 


1.8 ± 0.3 


1.0 ± 0.3 


7.9 ± 1.5 










cB58 


0.24^ 


1.7 ± 1.0^ 


12.6 ± 0.4^ 


1.1 ± 0.3^ 


< 2.3 


< 1.3 








8 O'clock 


0.57 




17.6 ± 0.5^ 


4.8 ± 0.4^ 


11.5 ± 2.2 











Note. — Col. (1) 24 /im fiux in mJy. /(24 //m) of Abell galaxies except Abell 2667a are from IRigbv et al.l (|20Q8I ). who 
noted that DAOPHOT errors of the fiux a re overl y optim istic and thus not reported. A nominal uncertainty of 0.1 mJy is 
adopted for SED fitting for L(TIR). Rigbv et al.l ([2008j) mistakingly reported a /(24 /im) for only one component of the 
A2667 arc. A correct value is reported here; Col. (2) 70 /im fiux in mJy; Col. (3)- (4) optical line fiuxes in 10"-*^^ erg/s/cm^; 
Col. (5)- (6) IR recombination lines fiuxes in 10 erg/s/cm^. The formal fit to the Bra line of cB58 yields a fiux value of 
(—0.4 ± 1.3) X 10"-*^^ erg/s/cm^ and hence a la limit is reported in the Table. Also for the Paa line, the formal fit yields 
(—0.2 ±2.3) X 10~^^ erg/s/cm^ and hence the la limit is reported. Likewise, a la limit is reported for the Paa line fiux of The 
Clone, where a formal fit yields (0.3 ± 1.3) x lO-^^; Col. (7)-(8) Line fiuxes of [Ne II] 12.8 /im and [Ne HI] 15.5 /im in lO-^^ 
erg/s/cm^ measure d by PAHFIT. Dots in dicat e no spectral coverage. These values are not corre cted for lensing m agnification. 
Flux references: (a) IHainline et"all (|2009l ): (b) ISiana et all (|2008l ): (c) ITeplitz et all (|2000l ): (d) IFinkelstein et al.l (2009). 
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Rest Wavelength (/xm) 



0.655 0.657 0.659 3.8 3.9 4.0 4.1 4.2 4.3 4 6 8 10 12 14 16 




1.332 1.336 1.340 7.8 8.0 8.2 8.4 8.6 10 15 20 25 30 35 



0.655 0.657 0.659 1.80 1.85 1.90 1.95 2 4 6 8 10 




2.303 2.310 2.317 6.4 6.6 6.8 10 15 20 25 30 35 



0.655 0.657 0.659 1.80 1.85 1.90 1.95 2 4 6 8 10 




2.335 2.342 2.349 6.4 6.6 6.8 7.0 10 15 20 25 30 35 

Observed Wavelength (/xnn) 



Fig. 3. Same as Figure [2] but for Abell 2667a, Abell 2218a, and Abell 1835a, from top to bottom, respectively. The IRS SL2 module 
spectra (not available for Abell 2218b in the previous figure), which cover Pao; lines for Abell 2218a and Abell 1835a, are shown in blue. 
There is no SLl and LL2 coverage for Abell 2218a and Abell 1835a. 



TABLE 4 
Measured PAH Fluxes 



Source 


/(PAH6.2) 


/(PAH7.7 Complex) 


/(PAHg.e) 


/(PAH11.3 Complex) 


/(PAH12.6) 


Abell 2218b 


15.9 ± 0.2 


72.2 ± 1.8 


13.2 ± 1.0 


17.3 ± 0.7 


12.7 ± 0.9 


Abell 2667a 


13.0 ± 0.2 


82.1 ± 3.0 


19.9 ± 2.2 


21.6 ± 1.3 


11.7 ± 1.6 


The Clone 


7.5 ± 0.3 


15.9 ± 0.4 


4.2 ± 0.3 


4.5 ± 0.7 




Abell 2218a 


16.0 ± 0.8 


54.4 ± 1.3 


10.3 ± 0.8 






Abell 1835a 


7.5 ± 1.4 


26.1 ± 1.6 


7.5 ± 1.9 






cB58 


2.9 ± 0.2 


9.7 ± 0.5 


2.5 ± 0.4 







Note. — PAH emission lines and emission complexes fiuxes measured by PAHFIT in 10"-*^^ 
erg/s/cm-^. These values are not corrected for lensing magnification. Dots indicate that the feature 
is outside spectral coverage. The 7.7 /im and 11.3 /im complexes are the sum of fiuxes at wavelength 
7.3 — 7.9 /im and 11.2 — 11.4 /im, respectively. The LL2 spectrum of the 8 O'clock arc does not 
have sufiicient S/N to measure PAH line fiuxes consistently. 
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Rest Wavelength (yam) 

1.8 1.9 2.0 2.1 2 4 6 8 10 




6.5 7.0 7.5 8.0 10 15 20 25 30 

Observed Wavelength (/xnn) 



Fig. 4. — Our reduction of the Spitzer IRS spectra for the 
Clone (SDSSJ1206+5142) and cB58 (cf. Fadelv et al. (2010) and 
[Siana et al. (2008), respectively). The color coding for IRS mod- 
ules the line fluxes for each Av scenario is the same as Figures [2] 
and [3] We note that the apparent spectral feature at ~ 5.9 /xm 
(observed frame) in the SL2 spectrum of the Clone cannot be in- 
terpreted as a Pao; line at the redshift z = 2.00, which has been 
measured from multiple optical lines (see discussion in Section [3?!]) . 

details in Section ISTTl The lensing m agnification of cB58, 
- 30x, was from Seit z et alJ (jl998f ). 

The Clone, at reds hift z = 2.0 026 was discovered from 
SDSS imaging (Lin e t al.l l20Q9l). Its rest-fr ame optical 
spectroscopy was obtained bv' Hainline et al.l (j2009). IRS 
ob servations in t he LL 2 and LLl bands were obtained 
bv iFadelv et al.l ()2010f ). who concluded that the Clone's 
IR emission is dominated by a starburst based on the 
strength of the 6.2 /im PAH feature, despite the object 
showing a very strong [SIV] emission line. The Clone's 
deep IRS SL2 observation is a part of Spitzer program 
50372 (PI C. Papovich) with an integration time of 8.6 ks. 
As with cB58, we re-reduced all IRS data using our pro- 
cedure and software to ensure consistency. We present 
our reduction of these archival data in Figure IH The lens- 
ing magnification for the Clone, 27 ± Ix, was reported 
bv lLin et all (120091 ). 

The 8 O'clock arc, at redshift z = 2.7308 was discov- 
ered s e rendipitously from . SPSS imaging bv lAllam et al.l 
(j2Q07[ ). iFinkelstein et al.l (j2009f ) obtained the rest-frame 
spectrosc opy of compo n ents ^ A2' and 'A3' in the nomen- 
clature of lAllam et al.l (j2007l ). The arc was observed as 
a part of Spitzer program 40443 (PI C. Papovich) with 
an SL2 integration time of 14.4 ks. The slit for the SL2 
observation is centered on the 'A2' component and cov- 
ered significant parts of both the 'Al' and 'A3' com- 
ponents such that any bias due to spatial variation of 
recombination lines along the lensed arc will be negligi- 



Rest Wavelength (yu.m) 




0.0 b I I I I i_ 

6.6 6.8 7.0 7.2 7.4 

Observed Wavelength (fim) 



Fig. 5. — The Spitzer IRS SL2 spectrum for the 8 O'clock arc, 
showing the Paa line. The color coding for IRS module and the 
line fluxes for each extinction scenario are the same as Figures [2] 
[31 and HI 

ble. The PAH spectrum of the 8 O'clock arc shows 6.2 
fim and 7.7 /im emission features, yet at low S/N due 
to contamination in the IRS slit from a nearby IR-bright 
dust-obscured galaxy, which we will discuss in the next 
section. We adopt the magnification estimate of 8x for 
the arc (|Finkelstein et al.ll2009f ). 

2.6. Additional Spitzer MIPS Photometry and Data 
Reduction 

The Spitzer MIPS 24 and 70 /im photometry for 
objects behind t he Abell clusters are taken from 
iRigbv et al] (|2Q08') with exception of Abeh 2667a, whose 
70 /im pho tometry is reported as an upper limit in 
lEigby et al." (^2008) but which was re-observed with 
longer integration time in Spitzer program ID 50586 (PI 
G. Rieke). For the Clone, we downloaded the 24 and 
70 /im data observed as part of programs 40430 (PI S. 
Allam) and 50372 (PI C. Papovich). The Spitzer MIFS 
observations for the 8 O'clock arc were obtained under 
program 40443 (PI C. Papovich). 

We measured the 24 and 70 /im fiuxes from the 
Level -2 (PBCD) data us in^ the Star finder IDL rou- 
tine (jPiolaiti et al.l |2QQQ[ ) to perform PSF photometry. 
The PSF model was generated using the STinyTim rou- 
tine provided by the S SC and smoothed according to a 
prescription given by Enge lbracht et al.l ((2007) to bet- 
ter match the observed PSF. We require a = 3 and the 
PSF correlation value of 0.75 for our 24 and 70 /im fiux 
measurements. The 8 O'clock arc has a 24 /im-bright 
dust-obscured galaxy ~ 6'.' 5 away that partially blends 
with the arc, but we were able to deblend the fiux by 
performing PSF photometry to subtract the intervening 
source first and repeat the photometry to measure the 
fiux of the arc. This contaminating galaxy is highly ex- 
tincted, such that it is undetected in the SDSS imaging 
and was discovered only in the 24 jam observation. Our 
fiux measurements along with fiuxes from the literature 
are reported in Table [3l 

2.7. Augmented Sample Summary 

Our augmented sample consists of the four galaxies we 
observed plus the three from the literature. In addition 
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to the measurements of Ha and of at least one of the IR 
recombination lines (Paa and/or Bra), all of them have 
good 24 fim photometry, far-IR and sub-mm photometry 
as well as IRS spectroscopy covering, at least, the 6 — 
8 /im wavelength region. The IRS coverage goes out 
to nearly 17 fim (rest-frame) in the two lowest redshift 
objects (Abeh 2667a and Abeh 2218b, z - 1.0). 



3. RESULTS 

In this section, we describe the near- and mid-IR spec- 
tra and the derived physical quantities, including metal- 
licity, L(TIR) and SFR from optical and IR tracers, and 
extinction measurements. 

3.1. Near- and Mid-IR Spectroscopy 

The final ID Ha spectra from our LBT/LUCIFERl 
observations are shown in Figures [2] and [3] along with 
the IRS spectra. The 6.2 /im PAH emission line, 7.7 jam 
complex, and 8.6 /im complex are clearly detected in all 
galaxies. 

The spectra from both slit positions (Section 12. 2p of 
Abell 2218b show highly asymmetric line profiles. The 
profile at the central position is double-peaked in both 
Ha and [Nil] while the central+0'.'9 spectrum shows a 
single-sided profile with stronger emission at the red side 
of the line. The combined line profile from both slit po- 
sitions is well fitted with a double Gaussian, whose de- 
convolved spectral width (corrected for instrumental res- 
olution) indicates a velocity dispersion of 163 ± 3 km/s. 
Traces of the arc northeast of the main component hinted 
in the ACS image are not detected at Ha. The IRS spec- 
trum for Abell 2218b shows a strong [Ne HI] 15.56 /im 
line as well as the [Ne II] line at 12.81 /im, apart from 
other PAH features. The Bra line is detected at 4.8<j 
with line flux of 4.3 ± 0.9 x 10~^^ erg/s/cm^. 

The spectra of Abell 2667a were extracted from two 
apertures in the slit. Both of these positions are very 
bright in Ha. Referring to Figure [U these apertures are 
at (1) the main arc north- northeastern of t he brightest 
cluste r galaxy (BCG) (identified as 'Al' bv I Sand et al.l 
(j2005'), and in Figure [1]); and (2) the smaller arc east 
of the BCG (the 'A2' arc in the Sand et al. study). Al- 
though it appears that a part of the smaller arc is outside 
the slit in the optical image, the Ha line is clearly de- 
tected (the light grey spectrum in Figure [3] top row, left 
column). Both apertures show the same profile for both 
Ha and [Nil] lines. No velocity shift between the two 
components is observed. The formal fit to the Bra line 
yields a flux of (0.4 ± 1.0) x 10~^^, a non-detection. The 
galaxy has the warmest dust continuum among our sam- 
ples, as seen in the rising SED towards the red end of 
the IRS spectrum in Figure [3] (cf. Abell 2218b's contin- 
uum in Figure [2]). However, the strong PAH emission is 
inconsistent with the system's emission being dominated 
by an AGN in the mid-IR. Its log([NII]/Ha) line ratio 
of -0.40 alone (i.e. without the [OIII] and H/3) does 
not unambiguously support or rule out the presence of 
a weak narrow-line AGN. But even with the presence 
of a weak, optically selected AGN, the EW of the 6.2 
micron PAH emission could still serve as a good tracer 
of SFR because the suppression of PAH emission due 



to AGN is found to be limite d to the nuclear region 
()Diamond-Stani"c fc Riekell20T(l . 

The Ha line of Abell 2218a also shows a double-peaked 
profile indicating a velocity dispersion of 188 ± 5 km/s. 
The galaxy has been a subject of a previo us near-IR spec- 
troscopic study by Richard et al] (|2011l ). Our slit orien- 
tation, however, is different (see their Figure 1); they 
excluded the /3 and 7 components (the for mer being a 
major component) in the nomenclature of iKneib et al.l 
()2QQ4l ). shown in Figure [U while our slit orientation 
covers all the main components (a, /3 and 7). The 
Paa observation, which targets the /3 component, de- 
tected Paa line at 4.8<j with line flux of 5.7 ± 1.2 x 
10~^^ erg/s/cm^. This measurement is consistent with 
the previously p ublish ed reduction of the same data by 
iPapovich et al.l (|2009f ) who measured the line flux to be 
8.5 dz 1.4 X 10"-^^ erg/s/cm^. We attribute the smah 
difference in the two mea surements to th e diffe rent def- 
inition of the continuum: IPapovich et al.l (|2009l ) defined 
their continuum with IRAC photometry whereas our re- 
duction adopted the simultaneous mul ti-co mponent fit of 
spectral features described in Section |2^ 

Abell 1835a shows a faint Ha line, too faint to mea- 
sure asymmetry in the line profile, but a bright Paa line 
with a flux of 7.9 ± 1.5 x 10"^^ erg/s/cm^ This is 
the second brightest Paa detection in our sample, after 
the 8 O'clock arc, whose Paa flux is 11.5 ± 2.2 x 10"^^ 
erg/s/cm^. The LLl PAH spectrum of the 8 O'clock 
arc does not have sufficient S/N to use PAHFIT to con- 
struct a noise free spectrum as a template for the Paa line 
flux measurement (Section l2.4p due to the contamination 
from the nearby IR-bright dust-obscured galaxy, thus in- 
stead the continuum of the arc for the measurement is 
defined by interpolating between the IRAC 5.0 and 8.0 
/im photometric points. These IRAC observations yield 
fluxes of 74 ± 3 and 64 ± 4 /iJy for the 5.8 and 8.0 /im 
bands, respectively. We present the Paa spectrum of the 
8 O'clock arc in Figure [51 

For the two IRS spectra from the literature, the Clone 
and cB58, our reduc tion generally agrees wit h those pub- 
l ished previously by iFadelv et al l (|2010[ ) and iSiana et al.l 
()2008f ). respectively. Small differences are noted in the 
PAH line fluxes due to the way that the continuum 
for the fit is defined. PAHFIT uses multi-component 
dust features to represent the continuum, while a single- 
component power-law is used by Fadely et al. and Siana 
et al. Their approach resulted in a similarly good fit to 
the data compared to ours, although the measured PAH 
flux could differ, as observed in the Clone' s PAH fluxes 
from t he 6.2 — 11.3 /im features where the iFadelv et al.l 
(201(f) values are a factor of 2.6 — 3.7 larger than ours. 
Fadely et al. reported tha t their results are a factor of 
1 — 8 higher than those of iBrandl et al.l ()2006[ ) in their 
re-reduction of the same data set. Brandl et al. use a 
continuum level deflned on either side of the emission 
features and thus their PAH fluxes will be systemat- 
ically smaller than with either our PAHFIT measure- 
ment or the single component continuum method. That 
is, we expect the measurements of the PAH line flux 
from the same data to yield flux values in the follow- 
ing increasing order: Brandl et al., PAHFIT, and Fadely 
et al. /Siana et al., which is consistent with the differ- 
ences we flnd. The formal Paa line flt for cB58 yields 
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TABLE 5 
Derived Quantities 



Source 


LIR (logLo) 






SFR (Mo/yr) 


* Ha/ Spitzer 
Av 




Av 




24 12m Far-IR SED 


Ra 


Ra + 24 12m 


PAH6.2 am LIR24 am LIRpar-lRSED 










(1) (2) 


(3) 


(4) 


(5) (6) (7) 


(8) 


(9) 


(10) 



Abell 2218b 


11.69 


11.68 


10.3 


36.6 


66.7 


48.0 


46.8 


2.8 ± 0.3 


2.1 


2.5 


Abell 2667a* 


11.68 




13.9 


39.9 


19.4 


47.1 




-2.4* ± 3.1 


1.7 


-0.4 


The Clonet 


11.20 


11.38 


11.4 


20.0 


36.2 


15.1 


23.2 


< -l.lt 


1.0 


1.0 


Abell 2218a 


11.81 


11.83 


7.5 


42.3 


170.2 


63.8 


66.5 


3.3 ± 0.4 


3.0 


3.2 


Abell 1835a 


12.71 


12.83 


14.5 


298 


558.2 


561 


752 


5.9 ± 0.4 


5.4 


5.7 


cB58t 


11.04 


11.04 


13.6 


19.8 


26.6 


10.4 


10.4 


< 0.7t 


-0.4 


-0.4 


8 O'clock 


12.17 




71.5 


154 




149 




2.8 ± 0.3 


1.0 


1.9 



Note. — Quantities in this table are corrected for lensing magnification. Col. (1) L(TIR) estimated from the monochromatic 24 /im 
indicator from Ruiopakarn et al. (2012); Col. (2) L(TIR) measured by fitting SED to far-IR and sub-mm photometry. The val ues for 
Abell 2218b and Abell 1835a are from Rigbv et al. (2008). The values for Abell 2218a, the Clone and cB58 are fro m Finkelstei n et al.l 
|2011), Fadelv et al. (2010) and Siana et al. (2008), respectively.; Col. (3) Ra SFR estimated by the Kennicutt! f[998) formula (no 
extinction correction); Col. (4) Extinction-corrected SFR estimated using the Ho; and rest-frame 24 /im luminosity formula provided 
by Kennicutt et al. (20091); Col. (5) Col. SFR estimated by the PAH 6.2 /im complex luminosity; (6) SFR estimated by 24 /im 
monochromatic indicator provided by Ruiopakarn et al. (2012); Col. (7) SFR estimated using the Far-IR L(TIR) in Col. (2) and the 
relationship given by Ruiopakarn et al. (2012); Col. (8) Screen extinction Av in mag estimated with Ho; and IR recombination lines 
assuming case-B recombination scenario. The value for cB58 uses Ho; and the Pao; line, the Bro; limit for cB58 indicates an extinction 
limit of Av < 1.8 mag.; Col. (9) A^^ in mag. See Section [3.41 for more details.; Col. (10) Final extinction estimate from averaging 
values in Col (8) and (9). For cB 58 and The Clone, the A^^ values are adopted. 
* IR recombination line detected at low significance, 
t IR recombination line not detected. 



a flux of (-0.2 ± 2.3) x 10"^^ erg/s/cm^ from which 
we determined the Paa flux limit to be < 2.3 x 10~^^ 
erg/s/cm^, which agrees with the non-detection reported 
bv lSiana et gd] (^008), < 6 x 10"^^ erg/s/cm^. The for- 
mal Bra line fit yields a line flux of (-0.4 ± 1.3) x 10"^^ 
erg/s/cm^, thus we report a limit of < 1.3 x 10~^^ 
erg/s/cm^. 

The formal fit to Paa line flux of the Clone yields 
(0.3 ±1.3) X 10~^^ erg/s/cm^, from which we report an 
upper limit similar to that of cB58 above. In the panel 
showing the Paa wavelength region for the Clone in Fig- 
ure m we note a spectral feature at 5.9 jam (1.96 /im 
rest-frame) that resembles a Paa line at z ~ 2.1. How- 
ever, the spectroscopic re dshift of 2.0026 was measured 
with multiple optical lines (iHainline et al.|[2QQ9l ) and thus 
precludes the possibility of the feature being interpreted 
as the Paa line. It is also unlikely to be an H2 emission 
line at 1.96 /im given the absence of other II2 lines nearby 
(2.03 and 2.12 /im). The source of the apparent feature 
therefore remains unknown. 

3.2. Metallicity 

We estimate metallicity with the N2 index 
(|Pettini fc Pagell 120041 ) using both linear and 3'^- 
order fits. The index is not affected by the optical 
extinction given the proximity of the lines. Both estima- 
tors yield consistent values, which agree well with those 
found by Erb et al. (2006) for massive galaxies at z > 2 
selected via UV luminosity. The range of the metallicity 
found in our sample is also in good agreement with the 
sample of 28 l ensed galaxies at 1.5 < z < 5 observed by 
iRichard et aLl (;2011i) . 8.00 > Z > 8.94, with an average 
oxygen abundance of Z = 8.55, similar to the average 
value from the linear N2 index in this work, Z = 8.56. 
That is, our sample has typical oxygen abundance for 
field and lensed galaxies at similar z. 



TABLE 6 

Metallicity Estimates from [NII]/Hq; 



Source 


Zppo4 


2pP04, S'^d order 




(1) 


(2) 


Abell 2218b 


8.55 ± 0.04 


8.53 ± 0.10 


Abell 2667a 


8.67 ± 0.05 


8.73 ± 0.10 


The Clone 


8.49 ± 0.02 


8.45 ± 0.07 


Abell 2218a 


8.55 ± 0.10 


8.52 ± 0.25 


Abell 1835a 


8.75 ± 0.15 


8.91 ± 0.33 


cB58 


8.31 ± 0.10 


8.26 ± 0.35 


8 O'clock 


8.58 ± 0.04 


8.57 ± 0.09 



Note. — Values are in 12+ log(0/H). Col. 
(1) Metallicity estimates from the N2 linear 
formula provided by Pettini & Pagel (|2004l ) : 
Col. (2) Metallicity e stimates from the N2 3^ ^- 
order estimator from lPettini Pagell (|2004I V 



3.3. Infrared Luminosity and Star Formation Rate 

We can measure the star- formation rate (SFR) from 
L(TIR), which in turn can be estimated by three trac- 
ers: (1) far-IR SED fitting; (2) the mo nochromatic 2 4 
/im L(TIR) estimator given by Ruiopak arn et al.l (|2Q12f ): 
(3) PAH emission line luminosity. The SFR can also be 
estimated from the Ha luminosity as well as the com- 
bination of Ha and 24 /im luminosity. We adopt the 
SFR derived from the SED-fitted L(TIR) as a fiducial 
for comparison. SFRs a re calc ulated from LfTIR ) by 
the relationship given by iRujopakarn et al.l ()2Q12f ) . All 
SFR estimates are converted from their respective IMF 
assumptions to the lKroupal (j2QQ2f ) IMF. The SFR results 
are presented in this section and their implications are 
discussed in Section HJ 

The far-IR SED fitting (e.g., 30 /im to ~1 mm) is 
the most straightforward among the methods to measure 
L(TIR). However, far-IR photometry is often affected by 
confusion noise at longer wavelengths in crowded fields, 
thus limiting the applicability of far-IR SED fitting in 
deep cosmological surveys. The far-IR and sub-mm pho- 
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tometry in our sample is aided by magnification of the 
gravitational lenses and thus we can probe fluxes below 
the typical confusion limit of the instrument. The SED 
fitting for the L(T IR) for Abeh 2218b and Abeh 1835a 
were carried out by Rigby et al. (2008) by combining the 
24 and 70 /im data with submillimeter observations, in- 
cluding at least 450 and 850 jam for all their galaxies 
and additionally 1.3 mm for Abell 1835a. The value for 
Abell 2218a is ni e asure d using Herschel observations by 
[FinkelsteinJEaD (|20TTh . The S ED-fitted L(T IR) for the 
Clone and c B58 w ere from the iFadelv et al.l (pOlO) and 
ISlana etld] (j2008) studies, respectively. Abell 2667a and 
the 8 O'clock arc have no L(TIR) estimates from far-IR 
SED fitting because the longest band we presently have 
is 70 /im. The 70 /im imaging for the latter is compli- 
cated the a nearby bright dust-obscured galaxy (Section 
12. 6p that blends with the arc. 

The monochr omatic (i.e. single - band) 24 /im L(TIR) 
estimator from iRujopakarn et all (j2012f ) yields L(TIR) 
consistent with the far-IR SED fitting at < z < 2.8 
within 0.02 dex, on average, with a s catter of 0.12 dex 
for individual star- forming galaxies. IRujopaka rn et al.l 
(|2012) take into account the evolution of the bolomet- 
ric correction from local galaxies out to z = 2.8 (the 
farthest redshift where the 24 jam band traces predom- 
inantly dust and PAH emissions) due to the structural 
differences of IR galaxies by using the IR surface bright- 
ness, Sl(tir)7 as an indicator of IR SED, which has 
been demonstrated by IRujopakarn et al.l (|201lll . Specif- 
ically, high-z IR galaxies are physically extended and 
typically have 100 — 1000 x more sur face area than lo- 
cal c ounterparts at similar L(TIR) (jRujopaka rn et aP 
l2011f ). and thus larger PAH-emitting area, which affects 
the bolometric corrections in the PAH wavelength re- 
gion th at is probed by the 24 / im band, particularly at 
z > 1 (jRujopakarn et al]l2012f ). We report the values 
of I/(TIR) from this 24 /im indicator along with those 
from SED fitting in Table [5l The agreement between the 
monochromatic indicator and the far-IR SED fitting is 
within 0.12 dex (average differ ence of 0.06 dex), simila r 
to the agreement reported in iRujopakarn et al.l (|2012l ) . 
No systematic trends in redshift, luminosity, or metallic- 
ity are observed. 

The last I/(TIR) estimator is derived via the luminosity 
of the 6.2 /im PAH emission line, following the relation- 
ship given by IPope et al.l ((2008|), which was measured 
by combining a local sample from iBrandl et al.l (|2006) 
and the res u lts for their sub- millimeter galaxies (SMGs). 
iPope et aH (|2008f ) report a slope of the relationship be- 
tween log 6.2 /im PAH luminosity and log L(TIR) to be 
approximately unity (formal slope = 1.01 ±0.01) over 
I/(TIR) ranging from 10^° - 10^^ Lq. Our objects are 
shown in Figure [6] on top of data points from iPope et all 
(I2008D and .Brandl et al.l (j2006f ). The 24 /im L(TIR) val- 



ues are used here instead of those from the far-IR SED fit- 
ting because Abell 2667a has no far-IR I/(TIR) estimate. 
Because of the methodolo fflcal differences in the PAH 
flux measurements between IPope et all ()2008l ) and ours, 
particularly the assumption of the continuum level, for 
the purpose of the comparison in Figure [6l we normalize 
the PAH luminosity by a factor determined by matching 
our 6.2 /im PAH luminosity of cB58 to the measurement 
of the same object by Siana et al.l (j2008f ). who adopt the 
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Fig. 6. — The relationship of 6.2 /xm PAH emission Une lumi- 
nosity to L(TIR) in our sample (large red squares) compared to 
the local star-forming galaxies and LIRGs from the IBrandl et al.l 
(1200 6") study (squares) and the high-redshift sub-mm galaxies from 
IPope et al.l (^2008) (crosses). The dashed line is a linear relation ship 
betwe en log luminosities with slope of 1.01 measured by Pope e t al.l 
("2008") . Our galaxies fall within this relationship, confirming the 
approximately unity slope reported by Pope et al. at the interme- 
diate luminosities between the local and SMG samples. 



same continuum assumption as iPope et all (^2008^. The 
6.2 /im PAH feature was chosen to represent the aro- 
matic emissions over other stronger PAH emission lines 
and emission complexes because it is narrow and rela- 
tively isolated from other lines, which mitigates the sys- 
tematic uncertainties introduced by the process of de- 
blending contributions from adjacent lines (which is re- 
quired for, e.g., the 7.7 /im complex and 8.6 /im line). 
To estimate the SFR from the 6.2 /im PAH luminosity, 
the PAH lu minosity is first converted to L(TIR) via the 
'Pope et al. l (l2QQ8i) conversion and then to SFR via the 
Rujopakar n et al.l (2012) relation. 

In addition to estimating SFR via L(TIR), we also 
show: the Kennicutt (1998) formula for Ha ( no extinc- 
tion correction) and the Ke nnicutt et al.l (j2009f ) indicator 
that combines Ha and the rest-frame 24 /im luminosity in 
order to correct for dust extinction in the SFR estimates. 
The SFR estimates from each method are tabulated in 
Table [5l We discuss the effects of extinction on the SFR 
estimates in Section HI 

3.4. Extinction 

Extinctions derived from optical lines (e.g. Ha/H/3) 
could suffer systematic underestimation where the over- 
all extinction is large and the dust is mixed with the 
sources. This possibility can be tested with IR recom- 
bination lines which are mo re robust to high extinction 
(|Alonso-He rrero et al."2006V To obtain a fiducial esti- 
mate, we assume a foreground screen and measure the 
extinction relative to Ha using the Paa and Bra lines 
assuming case-B recombination. For instance, for Bra, 



the A^*"/*^*^^'- 



is given by 



rQ!,0 



where /bf 



^,o/^Ha,0 



exp 



4 Ha/ Spitzer 



1.086 



-{ABra — ^Ha) 



(1) 



0.0291 and for Paa, Ipaa,o/lHafl 
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0.1226 (lOsterbrock fc Ferlandl 120061) for Te = 10,000 K 
at the low-density limit. The Aho^: Apan^, and A-^rcy Siie 
given by interpolation using the lRieke fc Lebofskvl([l985h 
extinction law, which for Av of 1 mag are 0.8 mag, 0.15 
mag, and 0.04 mag, respectively. 

Alternatively, the extinction can be estimated from the 
optical and IR SFR values by taking the latter to be a 
fiducial SFR. The assumption is secure in our sample of 
IR luminous galaxies, where direct UV leakage from the 
galaxy is expected to be small (Rieke et al. 2009 and 
references therein). This alternative extinction estimate, 
A^^, following IChoi et al.l (j2006h is 



aIR 



2 5 log f^E^ 

\SFRhq; 



(2) 



The A^ can then be converted to the Al^ via the 



iRieke fc LebofskyI (jl985f ) extinction law. When avail- 
able, we use the far-IR SFRjr to calculate the A^, oth- 
erwise the value estimated from the 24 /im L(TIR) is 
used. 

The Ar^^' method assumes a foreground screen 
of dust and therefore is a lower limit, whereas the A^ 
does not. The comparison of measurements from both 
methods, tabulated in Table [5l shows that in three out 
of four galaxies where IR recombination lines are well 
detected (Abell 2218b, Abeh 2218a, and Abeh 1835a), 
both methods agree within the range of uncertainty, 
which suggests that the nature of the dust distribution 
in these galaxies is roughly uniform, resembling the the 
foreground screen assumption. The other object with a 
well-detected Paa (the 8 O'clock arc), however, shows 
a 1.8 mag difference between a^^/^p'^^^^^ g^^id A^J^. The 
latter disagreement indicates an inhomogeneous mixture 
of dust in the 8 O'clock arc and highlights the diversity 
of the dust distribution scenarios at redshift 1 < z < 3. 

Abell 2218a, Abeh 1835a, cB58, the Clone, and the 
8 O'clock arc have extinction measurements from opti- 
cal and/or rest-frame optical spectroscopy in the litera- 
ture that can be co mpared with our IR measurements. 
iRichard et all (|2QTTI ) found £;(5-F)star of 0.18 for Abeh 
2218a from SED fitting, implying Ay of 0.6 mag, assum- 
ing the R = 3.1 law. Abell 1 835a has extinction mea- 
surements bv lNesvadba et all (|2007[ ) using IIa/H/3 ratio 
of E(B -V) = 1.3 - 1.6, implying Av of 4.0 - 5.0 mag. 
iTeplitz et all (|2000f ) measured E{B - V) of 0.27 for cB58 
from the line ratio, imp lying Av of 0.4 mag. For 

the Clone, iHainline etld] (|2009i ) measured E{B - V) us- 
ing the Hq;/H7 line ratio to be 0.28, which implies Av of 
0.9 mag. The extinction of the 8 O'clock arc were mea- 
sured bv lFink elstein et al. (2009) using weighted mean of 
to Ha, H/3, and H7 line ratios be Av = 1.17 ± 0.36 mag 
(using only IIa/H/3 yie lds E(B — V)^^^ of 0.97, irnplying 
Av of 1.3 mag); and bv lDessauges-Zavadskv et all (|2011i) 
to be E{B - l/jgas = 0.30 ± 0.04 mag, implying Av of 0.9 
mag. From Table O the extinction measurements from 
IR recombination lines for Abell 2218a, Abell 1835a, and 
the 8 O'clock arc, where we have secure IR line detec- 
tions are 3.3 ± 0.4 mag, 5.9 ± 0.4 mag, and 2.8 ± 0.3 mag, 
respectively. While our sample size is too small to draw 
a general conclusion, the comparison suggests that the 
optical measurements may sample systematically lower 
extinction regions, in agreement with local LIRGs and 
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Fig. 7. — Optic al extinction . Ay, from this work compared to 
the distribution of IChoi et al.l (|2006l ) from the Spitzer First-Look 
Survey (FLS) with mean redshift of 2; = 0.8. Our points are shown 
in filled squares, except for the cB58 point (empty square), where 
we adopt the A^^ value as an extinction estimate (see Section [374)) . 
The Ch oi et all (p006) points are A^^ measured by comparing the 
IR SFR, which is taken as a fiducial value, with the optical SFR 
measured from the Ha, H/3, and [Oil] emission lines. These are 
shown in red, blue, and grey circles, respectively. The intrinsic 
spread of points towards higher extinction of the Choi et al. sample 
may be larger if the line non-detections (upward arrows), which 
are likely due to large extinction, could have been measured. The 
intrinsic spread of extinction at high-z is far wider than previously 
known from optical-based measurements. 



ULIRGs fe.g. JAlonso-Herrero et aTI[2QQ6l ). 

The A^"/'^^'*^^^^ values for Abell 2667a and the Clone 
are negative because the formal fits to the Paa and Bra 
line yield values (or flux limits) lower than those expected 
for the extinction of mag (as shown in the simulated Av 
in Figures [3] and llj ; the line is undetected), albeit with 
a considerable uncertainty. For cB58, the line flux for 

both Paa and Bra are negative and thus the ^^'^^^^'^^^^'^ 
can not be formally calculated and we report a l-cr upper 
limit as cB58's ^^'^^^-^'^^^^^ . Given the general agreement 
between a^^/^p'^^^^'^ ^^Yid A^J^ in a majority of objects 
with secure measurements of Paa and Bra, we adopt an 
average value of a^^^^^'^^^^^ and A^J^ as a representative 
value of Av (except for the cB58 and the Clone for which 
we adopt the A^ value). 

We present our final Av estimates as a function of 
L(TIR) from the 24 //m indicator in Figure [71 Again, 
the 24 /im I/(TIR) is used here because Abell 2667a and 
the 8 O'clock arc lack far-IR I/(TIR) estimates. Our 
me asurements are co mpared with the A^ measurements 
bv IChoi et al.l (^2006), who study the extinctions of ob- 
jects selected in the near-IR and mid-IR from the Spitzer 
First Look Survey (FLS), which have a mean redshift of 
2; = 0.8 and luminosities ranging from the sub-LIRG to 
ULIRG range. Their extinction measurements are ob- 
tained by comparing the SFR estimates from the Ha, 
H/3, and [Oil] line fluxes with those from IR luminosities. 
That is, they estimated the extinction using the ratio of 
SFR from the individual optical indicators to the IR SFR 
estimates, e.g. SFRhc,/SFRir, SFRr^/SFRir, etc. The 
Choi et al. distribution of extinction values shows a cor- 
relation between A^ and I/(TIR) with a formal fit of 
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Fig. 8. — Comparison of star formation rate (SFR) estimated from the uncorrected Ho;, Ho; + rest-frame L(24 /im), L(TIR) from 24 
/im indicator, and L(TIR) from the 6.2 /im PAH emission Une, to the SFR estimated from far-IR SED fitted L(TIR), which is taken as 
a fiducial value, as a function of L(TIR) and Av: shown in the left and right panels, respectively. The Ho; SFR indicator systematically 
underestimates L(TIR) with increasing discrepancy from the fiducial at larger L(TIR), while the rest are consistent to the fiducial within 
the range of uncertainties. A trend line for the extinction effect on the Ha SFR indicator in objects with a screen-like extinction scenario 
(Section [4|) is shown as the blue dotted line. The points representing the 8 O'clock arc are circled in grey to indicate that L(TIR) of the 
arc is estimated from 24 /im because the SED-fitted L(TIR) is not available. The deviation of Ho; SFR of the arc from the trend line 
in the left panel, however, is not due to the L(TIR) measurement, but rather the differences of dust distribution scenarios between the 
inhomogeneous mixture in the arc and the homogeneous distribution in the rest of the sample (see Section |3.4|) 



A^J^ = 0.75 xlogL(TIR) -6.35, shown as the orange hne 
in Figure [71 The scatter is ^ 0.8 mag; however, this 
value does not include cases where [Oil] is undetected 
(shown as lower limits in Figure [7|). Allowing for these 
limits, the scatter is > 1 mag. This fit and th e large 
scatter are similar to the result of lAfonso et al.l (|2003), 
who estimate extinction by comparing Ra and [Oil] SFR 
estimates to those from the 1.4 GHz radio continuum at 
a sim ilar mean redshift of z = 0.8. A num ber of authors 
fe.g.. lBai et all 120071 : iKocevski et al.ll2011[ ) find that the 
variations in extinction are so large that [Oil] can be 
difficult to detect in some galaxies that have high SFRs 
indicated at 24 /im, although with high quality spectra 
there a re relatively few ca ses where [Oil] is completely 
absent (jWeiner et al.ll2007[ ). 

We found the spread of our extinction values to be 
at least as large. Abeh 2218a, Abeh 2218b, and the 8 
O'clock arc lie within the uncertainties of the Choi et al. 
relationship. Abell 2667a and cB58, with nearly zero mag 
extinction are less obscured for their L(TIR), whereas the 
extinction of Abell 1835a is far above the relationship. 
Selection effects may have reduced the scatter in many 
previous studies (e.g., the need to have an [Oil] line of 
measurable strength). Our results indicate that galaxies 
up to z ~ 3 show a very large range of Av, as is the case 
in the local Universe. 

4. EFFECTS OF EXTINCTION ON STAR 
FORMATION RATE INDICATORS 

In Figure [HI we compare the SFR estimates from 
Ra (uncorrected for dust extinction). Ha + 1/(24 /im), 
L(24 /im), and the 6.2 /im PAH emission with those of 
L(TIR) measured by far-IR SED fitting as a function of 
L(TIRj (Figure [8] left) and Av (Figure [8] right). Abeh 
2667a is excluded from this comparison due to its lack 
of the far-IR SED fitted L(TIR) and the fact that the 
curvature of the lensed arc caused a fraction of Ra flux 
to fall off the slit (the position of the slit is shown in 



the bottom left panel of Figure [TJ a part of the 'A2' arc 
is outside the slit). Without the correction for dust ex- 
tinction, the Ha SFR indicator is clearly affected at high 
L(TIR) and extinction. For the cases that extinction 
behavior resembles a uniform dust screen (every galaxy 
except the 8 O'clock arc), there is a tight trend of the 
Ha SFR deviation from the fiducial SED-fitted SFR as 
a function of the L(TIR) that is given by 

log[SFR(Ha)/SFR(LIR)] = 11.21 - 1.01 log(LIR) (3) 

Applying a single overall extinction to Ha will increase 
the estimates of SFRs and improve the agreement to the 
L(TIR) SFR, but wih not correct the trend with L(TIR). 
Further improvement in Ha-based SFRs require intro- 
ducing corrections for extinct ion as a functio n of stellar 
mass (e.g., Moustakas et all 120061: JWeiner et al . 2007), 
or of the SFR itself (e.g.[Buai_eLaD|2Q05|); se^Section [1 

The deviation of the Ha point of 8 O'clock arc in Fig- 
ure [8] (left) indicates that a larger amount of Ha flux is 
escaping from the arc given its L(TIR), which could be 
due to a different dust distribution relative to the star 
forming regions. The indicator is usually applied with a 
single nominal level of extinction; in this case, the line 
in Figure [8] (left) would be shifted upward but the slope 
would not be corrected. 

For Abeh 2218a, Abeh 1835a, and the 8 O'clock arc 
where Paa line is well-detected, we apply the extinc- 
tion correction to the Paa line luminosity and compare 
the ratio of 1/(24 /im)/I/(Paacor r) with the relationship 
that ' Alonso-Herrer o et al.l (|20Q6r ) found for local LIRGs. 
PaDovich et al. (200i) made this comparison for Abell 
2218a and found the L(24 /im)/L(Paacorr) ratio to be 
c:^ 0.5 dex lower than the local relation, indicating that 
Abell 2218a has lower L(24 jam) than local galaxies of 
comparable L(Paacorr), while agreeing with those of lo- 
cal individual HII regions. Papovich et al. interpreted 
this difference as an indication that Abell 2218a har- 
bors extended star-forming regions similar to a scaled- 
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up local HII regions, rather than nuclear starburst like 
local U/LIRGs. In our analysis, we have reproduced 
the L(24 /im)/L(Paacorr) ratio measured by Papovich 
et al. using our reduction technique for Abell 2218a, 
log[L(24 /im)/L(Paacorr)] = 2.2 ± 0.2, and also found 
that the ratio for Abell 1835a an d the 8 O^clock arc 
are 0.9 — 1.1 dex lower than the lAlonso-Herrero et al.l 
([2006 ) relationship . This is consistent wit h the result 
in Table [5] that the iRujopakarn et all ()2012l ) single-band 
24 /im L(TIR) indicator yields I/(TIR) values agreeing 
with the far-IR SED-fitted L(TIR), which implies that 
these galaxies have extended star formation (see Sec- 
tion 13.31 for details of the indicator). That is, the lo- 
cal relationship of LfPaacor r) and 1/(24 jam) reported by 
lAlonso-Herrero et al.l ()2006l ) could have limited applica- 
bility at high redshift because the Paa line will system- 
atically have lower extinction for a given L(24 /am) due 
to the extended structure of star formation. 

All the other indicators give consistent estimates, 
within the expected errors of ~ 0.2 dex. We foun d that 
the Ua + L(24 fim) indicator bv lKennicutt et al.l ([2009) 
tends toward a smaller value of SFR at large extinctions 
or large L(TIR), particularly above 10^^ L©, as observed 
in Abell 1835a and the 8 O'clock arc. The SFR estimates 
from the 6.2 /im PAH emission do not show a systematic 
trend with L(TIR), although the scatter is larger than 
with other indicators. The scatter for individual galaxies 
is ~ 0.2 dex, similar to the scatter of L(TIR) values we 
have found in Rujopakarn et all (|20 12) at redshifts where 
the 24 /im band probes the PAH emissions. Since these 
are high S/N measurements of an individual (i.e. rela- 
tively non-blen ded) PAH emissio n line, this result sug- 
gests that the IRujopakarn et all (2012) SFR indicator 
has succeeded in correcting for the SED evolution. 



5. CONCLUSIONS 

We observed four strongly gravitationally lensed star- 
forming galaxies at 1 < z < 3 with Spitzer/IRS and the 
LBT/LUCIFER to obtain the mid-IR and near-IR spec- 
troscopy. These observations are targeted to cover IR re- 



combination lines, including Ha in the near-IR and Paa 
or Bra in the mid-IR. We include another three galaxies 
from the literature with similar suites of observations, 
yielding a total sample of seven galaxies. Our sample 
spans the redshift range of 1.03 — 2.73 and the L(TIR) 
range of 1.3 x 10^^ to 7.0 x 10^^ Lq. 

The IR recombination line ratios are used to measure 
extinction that can probe deep into the highly obscured 
star-forming regions and thus provide an unbiased mea- 
sure of extinction under the foreground screen assump- 
tion. Independently, we estimate the extinction by com- 
paring the optical and IR SFRs, a method that does 
not make the foreground screen assumption. The re- 
sults from the two methods are consistent in three out of 
four galaxies with good IR recombination line flux mea- 
surements, suggesting that the dust extinction in these 
galaxies is consistent with a foreground screen (i.e. uni- 
form dust distribution). However, in the fourth case, the 
extinction estimates from two methods disagree by 1.8 
mag, indicating a deviation from the uniform dust screen 
assumption, which suggests an inhomogeneous dust mix- 
ture. The extinction range of our sample (assuming a 
foreground screen) is ~ 0.0 — 5.9 mag, which is a larger 
spread than previously known for intermediate and high 
redshift galaxies based on measurements with optical 
emission lines. These results suggest a large diversity in 
both the extinction levels and dust distribution scenarios 
at high redshift. 

We compare the performance of various SFR indica- 
tors over the extinction range and flnd that substantial 
extinction corrections are required for the Ha-based SFR 
indicator. The remaining indicators (1) combined Ha 
and L(24 /im); (2) L(24 /im); and (3) PAH (6.2 /im) ah 
give estimates consistent to within the expected uncer- 
tainties of ~ 0.2 dex. 

We thank Brian Siana and Eiichi Egami for helpful 
discussions. WR thanks Alexandra Pope and Philip 
Choi for data points in Figures [6] and [71 respectively, and 
acknowledges the support from the Thai Government 
Scholarship. This work is supported by contract 1255094 
from Caltech/JPL to the University of Arizona. 
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